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ABSTRACT 

The Blume - Emery - Griffiths model is investigated by use of the cluster vari- 
ation method in the pair approximation. We determine the regions of the phase 
space where reentrant phenomenon takes place. Two regions are found, depending 
on the sign of the reduced quadrupole - quadrupole coupling strength ^. For nega- 
tive ^ we find Para-Ferro-Para and Ferro-Para-Ferro-Para transition sequences; for 
positive ^, a Para_-Ferro-Para+ sequence. Order parameters, correlation functions 
and specific heat are given in some typical cases. By-products of this work are the 
equations for the critical and tricritical lines. 
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1. Introduction 

The Blume - Emery - Griffiths (BEG) model has originally been pro- 
posed to explain some features of the thermodynamic behavior of He^ - He^ 
mixtures'^] and has subsequently been used to describe the critical prop- 
erties of other physical systems like multicomponent fluids ^1 and magnetic 

systems PI. 

Since its introduction, it has been studied by many authors using a lot 
of different techniques: mean field approximation!^! , high temperature se- 
ries expansion^, Monte Carlo methods'^!, renormalization group^^"^], effec- 
tive field theories'^"-'^-'^!, cluster variation method (CVM)!-*^^]. In the last few- 
years, some authors'^"-*^^! pointed out the occurrence of a reentrance in the 
ferromagnetic-paramagnetic transition surface. In this paper, we use the 
cluster variation method'-*^^! (in the new formulation by An'-*^^! , which resorts 
to the Mobius inversion) in the pair approximation, to determine the region 
of the phase space of the model in which reentrance occurs. Two differ- 
ent regions are found, depending on the sign of the reduced quadrupole - 
quadrupole coupling strength ^. In the negative ^ region, the shape of the 
phase diagram leads us to characterize two kinds of reentrance: the first one 
is given by a first order Ferro-Para transition followed by Para-Ferro and 
Ferro-Para second order transitions; the second one is analogous except that 
it does not exhibit the first order transition. In the positive ^ region, a new 
Para_-Ferro-Para-|- reentrant phenomenon is noticed. In the most significant 
cases we give plots of order parameters, correlation functions and specific 
heat. Comparison of our results with those obtained with other techniques 
is made. By-products of this work are the equations for the critical and 
tricritical lines. 

The paper is organized as follows: in Sec. 2 we write the CVM free 
energy for the BEG model in pair approximation and the equations for the 
order parameters and the correlation functions; in Sec. 3 we determine the 
equations for the critical and tricritical points, that will be used in Sec. 
4 to determine the reentrance regions; in Sec. 4 we also give specific heat, 
order parameters and correlation functions for some typical cases of reentrant 
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phenomenon; in Sec. 5 some conclusions are drawn. 

2. CVM free energy and equations for the BEG model 

The BEG model has hamiltoniant^^ 

Hbeg = -JY. ^^^i - ^ E ^^i^l + DY^Sf-BY^Si, (1) 

where z,j are lattice site labels, Si is the 2;-component of a spin 1 operator 
at site z, Yli{ij) ^ ^'^^^ ^^e nearest neighbors (n.n.), and J > 
(ferromagnetic case) . In the pair approximation of the CVM (corresponding 
to the Bethe approximationt^^'^^J) one introduces the order parameters j/i = 
{Si) (( ) denotes thermal average) and y2 = (/Sf) and the n.n. correlation 
functions j/a = (SiSj), = {SiS^) and j/s = (SfSj) and considers the 
average value of the reduced hamiltonian 

n = iE^^jM. = dy, - 1(2/3 + M - byi, (2) 

zJ z 

where b = B/zJ, d = D/zJ and ^ = G/J are adimensional parameters and 
z is the coordination number of the lattice. 

The .educed foe energy pe. site / = = „ - ^| (where S is 

zJ zJ N 

the entropy) has been calculated in our approximation in Ref. 12, following 
An's scheme'-*^^!, in terms of the site and pair density matrices 

/ = 7i + r [(1 - z)Tr:{ps Inp,) + |lV(pplnpp)] , (3) 

with T = kBT/zJ. ps and Pp are respectively the site and pair density 
matrices, that turn out to be diagonal with elements 

y2 + yi 2/2- 2/1 ... 

Psi = — ^ — , = 1 - 2/2, Ps3 = — ^ — ; (4) 

_ ys + 22/4 + 2/5 _ 2/3 - 22/4 + 2/5 

PPI ~ ^ ' PP9 — ^ 5 

_ _ 2/2 - y5 + 2/1 - 2/4 _ _ 2/2-2/5-2/1+2/4 

Pp2 — Pp4 — 2 ' ^P*^ ~ ~ 2 ' 

Ppz = Ppr = = 1 + ys - 22/2- (5) 
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It is easy to verify that the relations 



Trps = Trpp = 1, ps = Tr pp 



(6) 



hold (the last Tr must be intended as a partial trace over the pair degrees 
of freedom, keeping fixed the site ones). In addition, the |/i's must obey the 
constraints given by 

0<Ps„Pp,<l (7) 

and have to be determined minimizing the free energy /. 

The minimum-/ requirement leads to the system of nonlinear equations 



df ^ 
— = 0, 

dyi 



i = 1, ...5. 



(8) 



This is a system in the yj's, that can be solved explicitly for the correlation 
functions but not for the order parameters. One can first solve the system 
for five of the Pp/s, obtaining 



' ppl 


= ViC^Pps 


PP2 




' PP3 


= iCVpp, 


PP6 


= U}2Vpp^ 


- PP9 


= mV'^Pp^ 



(9) 



where 
uji = exp 



d-b 



ZT 



UJ2 = exp 



r]i = exp 



l + ^-2{d-b) 



ZT 



C = 



y2 + yi 

[2(1-1/2) 



v = 



d + b 



ZT 



r]2 = exp 



2/2 - yi 



7 = exp 



-l + C-2d 



ZT 



l + ^-2{d + b) 



ZT 



and a = 



z - 1 



.2(1-1/2: 

Now can be determined by means of the condition Trpp = 1, which gives 
PP5 = W~^, with 



W = rjiC^ + r]2V^ + 2-fCV + 2uiC + 2^21^ + 1- 
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Finally, solving for the j/j's yields two coupled equations for the order param- 
eters Hi and j/2 

_ uJiC - UJ2V + ryiC^ - 772 

" (10) 
_uJiC ^uj2V ^r]iC'^ ^r]2V'^ + 2^CV ' ^ ^ 

and the expressions for the correlations functions 



2/3 = 



W 



VA = - — — • (11) 



w 



W 

It can be easily checked that eqs. (10) are substantially equivalent to the 
equations of Kikuchi's Natural Iteration Method^^^] , and thus they can be 
solved by iteration. Furthermore, they involve a reduced set of parameters, 
that is, the order parameters yi and 1/2, instead of the whole set of pair 
density matrix elements pp. , i = 1, . . . 9 contained in Kikuchi's equations. 

3. Stability and multicritical points 

In this section we derive the equations for the lines of critical and tricritical 
points, that will be useful in the next section to determine the regions of the 
phase space where reentrance occurs. We begin with the determination of 
the stability regions for the paramagnetic phase, that will lead us to the 
equations for the lines of critical points. 

Equations (10) and (11) are not sufficient to satisfy the /-minimum re- 
quirement, because they have been derived from the stationarity conditions 
(8). It is a well-known fact of analysis that / has a minimum in {yi} if 
and only if the Hessian matrix H = H({yi}) with elements Hij ~ — — - — 
is positive definite. In turn, H is positive definite if and only if det Hi > 0, 
i = 1, . . . 5, where 

Hii ■■■ Hii 



Hi = 



Hii ■ ■ ■ Hii 
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We have solved the inequalities detffj > 0, i = 1, . . .5 for the param- 
agnetic (that is, j/i = j/4 = 0) solution of (lO)-(ll) in absence of magnetic 
field (from this section on, we will take 6 = 0, and then uji = u>2 = uj, 
Vi = ^2 — v): ill order to determine the stability conditions for the param- 
agnetic phase. The first three inequalities are identically satisfied because of 
the constraints (7) on the density matrix elements. The fourth one gives the 
condition 1/2 > {z — 1)2/3. Solving the associated equation yields 



(12) 



which is satisfied by the second order critical temperature, where 

1 1 

C = e^/^"^ cosh — , 7o = (-S — 1) tanh — . 

ZT ZT 

Finally, the fifth inequality adds the condition 

{z-2)yl + y2-iz-l)y5>0, (13) 

which is identically satisfied for ^ < — 1. 

Plotting, for a fixed value ^ > —1, the solutions to eq. (12) and to the 
equation corresponding to (13), one obtains, if 2 > 3 (if 2; = 2 one has no 
phase transition), the diagram shown in Fig. 1. The phase space is thus 
naturally divided in five regions Si,i = 1,...5. Numerical analysis shows 
that in the region Si no stable paramagnetic phase can exist, in ^2, Ss and 
^4 there can exist only one stable paramagnetic phase, and in two stable 
paramagnetic phases can exist'^1, one with y2 < 1/2, that will be denoted 
by Para_, and another with j/2 > 1/2, that will be denoted by Para+. In 
this region there is a first-order Para_-Para+ transition surface ending at the 
line of critical points C^^\ at which y2 = 1/2. The behavior of the order 
parameter y2 around C is shown in Fig. 2. 

The equations of this critical line can be determined analytically by im- 
posing that the two solutions to the equation associated with (13) coincide. 
One thus obtains, choosing r, normalized temperature, as independent vari- 
able 



ln2 - ;2ln ( 1 - - 



z 



^^ZTln^ ^ . (14) 

cosh — 

ZT 
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Within these scheme it is also possible to determine the equations for the 
tricritical line, which separates second order transitions from first order ones. 
Indeed eqs. (10) can be written, to the 2nd order in yi, in the form (the 
solution j/i = has been factored away) 



Pq + P2yi = 

qo + q2yl = o 



(15) 



where 



a 



po = po{y2) = 2{ri + '^)Vo^ + AcoVq + l-2Vo — {io + 2riVo), 

2/2 

qo = 90(2/2) = 2(r/ + 7)2/2^0^ + 4a;2/2Vb + 2/2 - 2{r] + -f)V^ - 2a; Vb, 



P2 =P2{y2) = 



= 2 



2a 



^Voa + 2uj 



a 



2/2 



a\ ^/2a 
uj\ „ I + r/Vo 



q2 = 92(2/2) = 



2/2 



with Vb = X 
[2(1-2/2) 
and (15) becomes 



2a 
2 

2/2 



7VbQ; + 2uj 



a 



2^ 
2/i 



a \ (2a 



7VbQ: 



. Again, to the 2nd order in yi, we have j/2 = a + byf 



po{a) + {p2{a) + bp'o{a))yl = 
qo{a) + (q^ia) + bq'o{a))yf = 



(16) 



Letting yi = in (16) gives the second order critical temperature again, while, 
requiring that (16) are identically satisfied, we have the following system of 
equations for the tricritical line^^^^ 




(17) 



with all functions evaluated at a (which, of course, must be determined from 
one of the equations). Eqs. (17) cannot be solved analytically for two of the 
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three parameters ^, d, r, but one can solve them numerically at fixed ^ by 
looking on the second order critical line given by eq. (12) for the point at 
which the last of eqs. (17) changes sign. The values of r and d so obtained are 
plotted in Fig. 3a-b for typical values of the coordination number. Our results 
are compared, for 2; = 6, with the effective field theory ones by Tuckerl-*^^] in 
Table I, showing that our temperatures are slightly lower, with an agreement 
to within 5%. 

Finally, let's notice that in the limit r — > the critical and tricritical lines 
converge in ^ = — 1, (i = 0. 

4. Reentrant phenomenon 

At fixed — 1 < ^ < 0, our calculations show that the phase diagram has 
the general shape given in Fig. 4. The solid line represents a second order 
transition, while the dotted one represents a first order transition (determined 
numerically, by comparison of free energies). We have denoted by dt, do 
and dm respectively the value of d at the tricritical point, at r = 0, and 
the maximum value of d for which a transition can occur. It is evident 
from this diagram that we can have two different kind of reentrances: when 
do < d < dm we have, increasing the temperature, a Para-Ferro transition 
followed by a Ferro-Para one (both second order); when dt < d < min((io, dm) 
we have a Ferro-Para-Ferro-Para sequence of transitions, the first of which 
is first order, in contrast with the result of Ref. 18, and the others second 
order. 

We have already mentioned how dt can be determined in our scheme, 
while dm can be easily evaluated numerically by means of eq. (12). In order 
to determine do we must minimize the zero temperature free energy 

/o = ^ = %2- ^(2/3 + ^2/5). (18) 

This is a linear function of {yiiVziVb) defined in a polyhedron, and thus its 
minima must be located on the vertices of the polyhedron, which, according to 
the constraints (7), are Pq = (0, 0, 0), P1/2 = (1/2, 0, 0) and Pi± = (1, ±1, 1). 
The stable phase (that is, the one with lowest free energy) can be paramag- 
netic only if 1/2 = because if 2/2 > values of yi different from zero would 
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be allowed and a little magnetic field would make the phase with yi = Q 
unstable. This means that the stable phase is paramagnetic if and only if the 
absolute minimum of /o is located at Pq, that is, if d > do = 1" ^ ■ 

The values of do, dt and dm for 2; = 6 are plotted in Fig. 5, where the 
region of the d) plane where the two kind of reentrances occur are also 
indicated. For 2; < 4 we have no reentrance in agreement with Ref. 8 and 
in contrast with the results of Ref. 11, while for increasing z the reentrance 
regions become larger, as already noticed in Ref. 10. The maximum ^ for 
which there exist Para-Ferro-Para reentrance are ^ = —0.345, —0.283, —0.242 
for 2 = 6, 8, 12 respectively. We also give, for the Ferro-Para-Ferro-Para case, 
the order parameters, the correlation functions (Fig. 6a-b) and the specific 
heat (Fig. 7). 

If ^ < —1 there is no longer first order transition, because the tricritical 
temperature approaches zero as ^ tends to -1 (Fig. 3a). Again for < d < 
dm we have a Para-Ferro-Para sequence of transitions, but now do has to be 
determined in a different way: indeed, the transition near r = is second 
order, and thus we have 

do= lim 4(T) = l+e, (19) 

T-»0,^<-l 

where dcir) is the second order critical value of d, obtained from (12). It 
must also be noticed that in this case the low temperature paramagnetic 
phase might be a staggered quadrupolar phaset^^l. We do not give more 
details about this case, because it would require to consider two different 
sublattices, and thus to change the whole scheme. 

If < ^ < 3 we have no reentrance because both the second order and 
the first order transition lines have negative slope in the (d, r) plane. 

Let us examine the last case of reentrance. A numerical analysis shows 
that, for fixed C > 3, the first order Para-Ferro transition line in the (d, r) 
plane starts from r = 0, d = do = — ^-^ with positive slope, while the second 
order one has negative slope, as reported in Fig. 8 (phase diagram for 2; = 4, 
^ = 3.3). The point P is a four phase coexistence pointful, meeting point 
of three first-order transition lines, Para_-Ferro, Para-|--Ferro and Para_- 
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Para.). , and the corresponding value dp (that must be evaluated numerically) 
is the maximum value of d for which a Para-Ferro transition can occur. For 
do < d < dp we have a Para_-Ferro-Para+ sequence of transitions, and thus 
another reentrance. The values of do and dp are plotted in Fig. 9. We have 
checked that, even in this case, increasing the coordination number of the 
lattice the reentrance region becomes larger. Again we give plots of order 
parameters, correlation functions (Fig. lOa-b) and specific heat (Fig. 11). 

It is noteworthy that already the early mean field phase diagrams by 
Blume, Emery and Griflaths^^J could suggest the existence of this kind of 
reentrance. 

Finally we remark that the four phase coexistence point P and the tri- 
critical point T (Fig. 8), which are well distinct in mean field approximation 
and coincide in renormalization group studies!^], are very close in our ap- 
proximation. In fact, their locations in the ((i, r) plane are P = (2.17,0.75), 
T = (2.10,0.79) in mean field approximation, and P = (2.16,0.59), T = 
(2.15,0.60) in our approximation. 

5. Conclusions 

We have studied the Blume-Emery-Griffiths model in the pair approxima- 
tion of the cluster variation method, determining the phase space conditions 
for the reentrant phenomenon. Our results show that there are three differ- 
ent kind of reentrances, characterized by different transition sequences. For 
negative values of the reduced quadrupole-quadrupole interaction strength 
^ we have found that reentrance is allowed only for 2; > 6 and for ^ < ^m, 
where depends on z. For positive values of ^ we have found a little region 
of the phase space where a new Para_-Ferro-Para+ reentrance can occur. 
Both regions become larger for increasing z, as already pointed out in the 
literature. 

As by-products, we have given analytical expressions for the critical and 
tricritical lines. 
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Table I. Coordinates of the tricritical points in the {d,ZT) plane {z = 6). In 
the third column are reported results of Ref. 17 for comparison. 



d 


ZT 


ZTeft 


0.468 


1.61 


1.67 


0.506 


1.76 


1.85 


0.543 


1.92 


2.01 


0.580 


2.06 


2.17 


0.617 


2.19 


2.31 


0.654 


2.33 


2.45 
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Figure captions 



Fig. 1 : Stability regions for the paramagnetic phases at 2; = 6, ^ = 3.3. 
Fig. 2 : Quadrupolar order parameter j/2 around the critical point C (z = 
6, e = 3,3). 

Fig. 3a : Values of r vs. ^ at the tricritical point for 2; = 3, 6, 12 (solid lines) 

and in mean field approximation (dashed line). 
Fig. 3b : Values of d vs. ^ at the tricritical point for z = 6 (solid line, 

changing z the line does not vary significantly) and in mean field 

approximation (dashed line). 
Fig. 4 : Phase diagram for z ~ 6, ^ = —0.9. T is the tricritical point. 

The solid and dashed lines represent respectively second and first 

order transitions. 

Fig. 5 : Reentrance region at — 1 < ^ < 0. Solid line = do, dashed line 
= dm, dot-dashed line = dt- In the region P the ground state is 
paramagnetic and one has a PFP sequence of transitions; in the 
region F one has ferromagnetic ground state and FPFP transition 
sequence. 

Fig. 6a : Dipolar order parameter yi (solid line) and correlation function 
j/4 (dashed line), vs. r, for the FPFP reentrance {z = 6, ^ = —0.5, 
d = 0.2475). 

Fig. 6b : Quadrupolar order parameter 1/2 (solid line) and correlation func- 
tions (dashed line) and 2/5 (dot-dashed line), vs. r, for the case 
of Fig. 6a. 

Fig. 7 : Specific heat vs. r for the case of Fig. 6a. 

Fig. 8 : Phase diagram for z — 4, ^ — 3.3. T — tricritical point, C = criti- 
cal point, P = four phase coexistence point. The solid and dashed 
lines represent respectively second and first order transitions. 

Fig. 9 : Reentrance region at ^ > 3. Solid line = do, dashed line = dp. 

Fig. 10a : Dipolar order parameter yi (solid line) and correlation function j/4 
(dashed line), vs. r, for the P_FP+ reentrance (^ = 3.3, d = 2.16). 

Fig. 10b : Quadrupolar order parameter 7/2 (solid line) and correlation func- 
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tions Us (dashed line) and j/5 (dot-dashed hne), vs. r, for the case 
of Fig. 10a. 

Fig. 11 : Specific heat vs. r for the case of Fig. 10a. 
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